Time-resolved charge detection with cross-correlation techniques 
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We present time-resolved charge-sensing measurements on a GaAs double quantum dot with two 
proximal quantum point-contact (QPC) detectors. The QPC currents are analyzed with cross- 
correlation techniques, which enable us to measure dot charging and discharging rates for signif- 
icantly smaller signal-to-noise ratios than required for charge detection with a single QPC. This 
allows us to reduce the current level in the detector and therefore the invasiveness of the detection 
process and may help to increase the available measurement bandwidth in noise-limited setups. 
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The use of quantum point contacts (QPCs) as charge 
sensors integrated in semiconductor quantum dot (QD) 
structures"'^ has become a well-established experimental 
technique in current nanoelectronics research. The time- 
resolved operation of such sensorsSi^i^ allows us to ob- 
serve the charge and spin dynamics of single electrons^ 
which has potential applications in metrology^i or for 
the implementation of qubit readout schemes in quan- 
tum information processing.^ Another appealing prop- 
erty of the QD-QPC system is that it opens the possibil- 
ity of studying a well-defined quantum mechanical mea- 
surement process and testing the theory of measurement- 
induced decoherence^ 

The difficulty in achieving quantum-limited charge de- 
tection is mainly the limited bandwidth of the readout 
circuit compared to charge coherence times. In addi- 
tion, decoherence mechanisms exist that are due to the 
QPC but not directly linked to detection, such as the 
excitation of electrons in the QD driven by noise in the 
QPCfiS, an effect which is more pronounced at higher 
source-drain voltages. Both problems are related to the 
limit in signal-to-noise ratio (SNR) offered by present- 
day setups. A common experimental approach to over- 
come such a limitation is the use of cross correlation 
of independent measurement channels. In the context 
of charge sensing, correlation techniques have previously 
been used in Al single electron transistor setups to sup- 
press background charge noise^^ and to obtain estimates 
for the spatial distribution of sources thereof. High- 
frequency noise measurements usually rely on correla- 
tion techniques which eliminates noise contributions of 
the wiring and the amplifiers. ^"^ 

In the present work we present cross-correlated charge 
sensing measurements in a double quantum dot (DQD) 
sample with two charge readout QPCs. The potential 
advantages of such a design for the continuous quan- 
tum measurement of charge qubit oscillations have been 
put forward by Jordan and Biittiker,^* While the cor- 
responding time scales are yet beyond our experimen- 
tally achievable bandwidth, we demonstrate the bene- 



fit of cross-correlation techniques in the classical detec- 
tion of electron tunneling. By a detailed analysis of the 
cross-correlation function of the QPC currents and of 
higher-order correlators, we are able to measure tunnel- 
ing rates in a manner eliminating uncorrelated amplifier 
noise. Compared to a measurement of the same quan- 
tities using only one channel, we are able to reduce the 
detector current by roughly 1 order of magnitude. 
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FIG. 1: (Color onhne) (a) Inset: AFM micrograph of the 
sample which consists of two QDs in series (QDl and QD2) 
with two charge readout QPCs, denoted PCI and PC2. The 
source, drain, and center barriers can be tuned with in-plane 
gates S, D, and C. Main graph: part of the DQD charge sta- 
bility diagram obtained by counting the number of switching 
events in 72. (b) Detector currents recorded at two different 
gate configurations indicated in (a). Dot-lead tunneling pro- 
cesses (I) cause identical switching directions in both QPCs 
whereas inter-dot processes (II) cause opposite switching, 
(c) Current correlator Co = {hh) extracted 
from the raw data used in (a), revealing the correlation- 
ant icorrelat ion pattern in the Vs-Vd plane. Numbers (n, m) 
indicate the electron occupancy of the dots relative to the 
state (0,0). 
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The inset of Fig. [Ua) shows the structure, fabricated 
on a GaAs/AlGaAs heterostructure containing a two- 
dimensional electron gas 34 nm below the surface (den- 
sity: 5 X 10^5 and mobility: 40m2/Vs at 4.2 K). 
The electron gas was locally depleted by anodic oxi- 
dation with an atomic force microscope (AFM)?^^ The 
measurements were performed in a '^He/^He dilution re- 
frigerator with an electron temperature of about 80 mK, 
as determined from the width of thermally broadened 
Coulomb blockade resonances. ^-^ The structure consists 
of two QDs in series (denoted QDl and QD2) with two 
charge-readout QPCs (PCI and PC2). The strength of 
the tunneling coupling to source and drain leads is tuned 
with the gates denoted S and D; gate C controls the in- 
terdot coupling and is kept at a constant voltage for these 
measurements. 

Both QPCs are voltage biased and tuned to conduc- 
tances below 2e^//i. Their currents are measured with 
an I /V converter with a bandwidth of 19 kHz and sam- 
pled at a rate of 50kS/s. The data is stored for further 
processing in the form of time traces typically few sec- 
onds long. Electrons entering or leaving either dot cause 
steps in the currents that can be counted.^ Figure [TI^ a) 
shows a color plot of the count rate in PC2 vs S and 
D gate voltages close to a pair of triple points of the 
DQD system^^ at zero source-drain voltage. Lines with 
negative slope belong to equilibrium tunneling events be- 
tween the dots and the leads. The inter-dot charging en- 
ergy (0.3 meV) is much larger than the thermal energy, 
therefore also the line of inter-dot tunneling events with 
positive slope is observable. The corresponding tunnel- 
ing rate of about 1 kHz is the largest in the system. Few 
additional counts outside the main resonances are due to 
excitation processes driven by the currents in the QPCs 
(Ref. [13) (source-drain voltage 300 /iV) . 

Due to geometry, the capacitive coupling between the 
QPCs and the QDs is asymmetric; charging QDl will, for 
example, cause a larger step in the conductance of PCI 
than charging QD2. Accordingly, the steps due to dot- 
lead tunneling events have the same sign in both QPCs 
whereas inter-dot events cause opposite switching as seen 
in the time traces plotted in Fig.[TJb). A simple param- 
eter which characterizes the correlation between the two 
channels is the correlator 

^ {hh)-{h){h) 

° ^{r()-{h)N{i'i)-{h?' 

where angular brackets denote time averaging. We ob- 
tain this quantity, as well as any other cross-correlation 
expression discussed later in this paper, by digital pro- 
cessing of the raw time trace data. In Fig. [Ijc), we plot 
Co calculated from the same data as used in panel (a). 
It clearly displays the expected pattern of positive and 
negative correlations along the charging lines of the DQD 
stability diagram. Note that in the following, we implic- 
itly assume the mean values of Ii and I2 to be subtracted 
by setting (/i) = (/a) = 0. 

Going beyond this more qualitative information, in the 



following we analyze how to extract physical tunneling 
rates with the help of cross-correlation techniques and 
apply this to the example of tunneling from the lead into 
and out of QD2 (rates Fjn and Tout) in the present sam- 
ple. The underlying problem is to extract these two char- 
acteristic parameters of a random telegraph signal (RTS) 
which is, as we assume, a component of both QPC 
currents, along with uncorrelated noise. If the noise is 
stronger than the signal, the information on the actual 
time dependence of /'^'^^ is lost even if there are two mea- 
surement channels available. This is however not a prob- 
lem since one can determine the rates Fin, out entirely on 
the basis of time-averaged quantities derived from 
For the analysis presented here, these are on the one hand 
its autocorrelation function from which we can extract a 
characteristic time constant tq — l/(rin -|- Pout) and on 
the other hand its skewness 7 which depends on the oc- 
cupation probabilities of the high and low current states 
of l''^-' and allows one to determine the ratio Pm/Fout- 
The sought-after Pin, out are then uniquely determined by 
To and 7. This concept of exploiting third-order cumu- 
lants of a telegra ph s ignal for measurement has also been 
discussed in Ref. [H. 

To state this more precisely, we split up the QPC cur- 
rents according to Ij — ajl'^'^^ + j — 1, 2, where 
aj are dimensionless factors (ai > by convention) and 

/j"'' are mutually uncorrelated noise components. The 
product of Ii and I2 appearing in the cross-correlation 
function C(t) = {Ii{t)l2{t + r)) then consists of four 
terms among which any one containing a factor /^"^ or 
/j"'' is integrated to zero. The only nonvanishing part is 
then proportional to the autocorrelation function of the 
signal /(^^ , 

C(r) « aia2{l'^^Ht)I^"Ht + r)) 

= aia2(/(^)')e^l^l/^«, (2) 

where the decay time of the exponential is given by 
To = l/(Pin-Krout)'^ Note that form (O of C(t) implies a 
purely Poissonian tunneling process. On time scales rele- 
vant for our measurements, non-Poissonian statistics can 
occur when excited dot states are involved^*' and would 
manifest itself in a deviation of C (r) from the exponential 
shape. Figure mja) shows a set of C(t) curves belonging 
to the crossover from the (1,0) to the (1,1) state in the 
DQD charge stability diagram. For curves in the cen- 
ter of this plot, the electrochemical potential of QD2 is 
roughly aligned with that of the lead, and the tunnel- 
ing in and out rates are similar. The peak amplitude of 
C(r) is largest in this regime. It is proportional to (j''^)^) 
which is maximum in the case of a symmetric RTS, as 
we discuss later in more detail. Moving away from this 
point, the peak amplitude decays. The behavior of the 
peak width outside the resonance is determined by the 
behavior of the rates Fjn^out: While one of the rates tends 
to zero, the other approaches its finite saturation value 
which is also the saturation value of I/tq. The peak 
width therefore remains nonzero. 
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FIG. 2: (Color online) (a) Evolution of the current cross- 
correlation function C(r) = {Ii{t + r)/2(t)) when crossing 
the (1, 0) — > (1, 1) charging line in the DQD stability diagram 
[outside the scan range of Fig. mb)]. The curves are offset 
for clarity; going from the lowest to the uppermost curve cor- 
responds to adding one electron to QD2. C(t) exhibits an 
exponential decay ~ exp( — jrl/ro) characteristic for the ran- 
dom telegraph signal present in both QPC currents. The cen- 
ter curves, which feature the largest amplitude C(0), belong 
to QPC signals with rather balanced occupation of the high- 
and low-current states, (b) Noise reduction due to cross cor- 
relation in time traces measured on a Coulomb peak. Plotted 
are the Fourier transform Sc of the cross-correlation function 
and the geometric mean of the two QPC current spectra Spci 
and S'pc2- The dashed line represents the ideal RTS spec- 
trum 5'rts(/) = 2(7^r],g)ro/[l -I- tq (27r/)^] with parameters 
(-^RTs) = 0.22 X 10"^ nA^ and to = 0.8 ms. (c) Standard de- 
viation adhh)) of the average (/1/2) = T"^ h{t)h{t)At 
as a function of the integration time T. The solid line marks 
the expected ~ behavior. 



The noise reduction due to the cross correlation is 
best visualized in the frequency domain. In Fig. [D^b), 
we plot the geometric mean of the power spectral densi- 
ties of some example time traces Ii and I2 along with 
the Fourier transform of their cross-correlation func- 
tion. The spectrum of the raw traces consists of the 
Lorentzian contribution of the telegraph signal and a 
noise background on the order of 10~^nA^/Hz which 
is dominated by the (current-independent) noise of the 
room-temperature I /V converter and contains an addi- 
tional current-dependent part that is most likely related 
to charge noise in the sample. In the cross-correlation 
spectrum, the signal part is unchanged; the noise on the 
other hand is clearly suppressed. This remains, for the 
moment, a qualitative statement, and we postpone the 
quantitative discussion about the noise reduction to the 
end of the paper. 

The correlation time tq gives the sum of the two tun- 
neling rates but is insensitive to their relative magnitude. 
A second experimental parameter is therefore needed 
which depends on Fin/Fout and is also accessible in high- 
noise conditions. It is natural to consider the skewness 
7 = (/('=)3^/(/(c)2^3/2 i^gcause it parametrizes the degree 



of asymmetry in the current distribution function of /^"^^ , 
which is in turn fully determined by Fjn/Fout- Namely, 
the occupation probability of the low-current state of the 
RTS (electron on the dot) ispiow — rin/(rin+Fout); anal- 
ogously phigh = rout/(rin + Tout)- Assuming a current 
difference of A/ between the two states and (/'•^■') — 0, 
then the nth central moment of is given by 



= Plow(-PhighA/)" + PhighblowA/)" 
rinTout 



(rin + Tout) 



[F" 



]A/". 



(3) 



In calculating the skewness based on Eq. ^ for n ~ 2 
and 3, we see that the current scale A/, i.e., the informa- 
tion on the strength of the QD-QPC coupling, is elimi- 
nated. After some algebra, we obtain the expression 



7 ^ 



c)3\ 



(j(c)2)3/2 (Fi„Fout)l/2- 



(4) 



Using Eq. ([4]) and the previously determined tq = 
l/(rin + Tout); we can now write down the total event 
rate, 



I^in I^n 



Tin Tout 
(Ein + Fout) J 
1 



EinE, 



inJ- out 



^rinTout ~l~ Tq^^ 



7-0(4 + 72)' 

The individual tunneling rates are then given by 

2 



Ei 



in/out 



To(4 + 72 =F7v/4 + 7^) 



(5) 



(6) 



The skewness is experimentally accessed through an 
appropriate combination of second- and third-order cor- 
relators computed from the raw time traces Ii and I2 
that have the property to be insensitive to the back- 
ground noise. On the one hand, we use again (/1/2) 
which is the cross-correlation function at zero time differ- 
ence C(t = 0) and is equal to aiUj (J(=)2^. On the other 
hand, we use the combinations (iflj) which are propor- 
tional to the third moment of /(^^. In writing the QPC 



(") 



currents as a sum of signal and noise, , 

(n) 

it is readily seen that any term containing the noise Ij 
gives zero contribution to the time average and we have 
(iflj) ~ Q^fctji-^^^^'^)- The skewness can then be ex- 
pressed as 



7«sgn((/i/|)) 



1/2 



(7) 



The asymmetry in this formula is caused by our previous 
choice «! > 0; i.e., we fixed the sign of I^'^) such that it 
is positively correlated with Ii. This freedom of choice 
is not unique to our correlation analysis. Instead, the 
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FIG. 3: (Color online) Plots of the rates Fin, Tout, and 
Ftot — FinFout /(Fin+Fout) as determined by electron counting 
(marked with the letter "C") and by current cross-correlation 
( "X" ) . The sweep range is indicated by a line in the top graph 
of Fig.lHa). 



assignment of one detector event type (e.g., "PCI current 
up") to one system event type (e.g., "electron tunneling 
from QD2 into lead" ) has to be done in any case. 

Before turning to the experimental results, we discuss 
the role of the integration time T in the cross-correlation 
process. How large do we have to choose T until the 
cross-correlation function ([2]) is reproduced to a good ac- 
curacy? In order to estimate the remaining noise contri- 
bution to C(r) after averaging, we treat the integration 
as a summation over samples that are separated in time 
by the typical autocorrelation time of the noise r„ and 
are therefore statistically independent. Using the cen- 
tral limit theorem, we write the standard deviation of 

this sum as (T„/r)i/2(aJ"V^"V^' [cf- Fig.^a), inset], 
where we have introduced the symbols o-j"-* = ^j(")2^i/2 
for the noise in the channels. It should not exceed the 
contribution of the telegraph signal aia2 {l'-"^^). For the 
data presented here, time traces were recorded for 5 s 
and digitally low-pass filtered at 3 kHz (t„ « 0.1ms), 
yielding an expected noise reduction of 0.005. As seen 
from Eq. ([3]), the quantity aiQ;2(/*-'^''^) contains a factor 
rinTout and is therefore small whenever one of the tun- 
neling rates is small. As a result, the situation where 
the two rates are similar presents the optimal case for a 
correlation measurement. 

Even disregarding any uncorrelated noise, the expo- 
nential shape of the autocorrelation function of 7^'^-' is the 
limit of infinite integration time. It is practically reached 
under the condition that T covers a sufficient number of 
switching events, T ^ l/Ftot- This second condition on 
T is therefore linked to the statistical uncertainty of the 
measurement. 

In Fig. [3l we compare the outcome of the conventional 
(counting) method and the correlation procedure for a 
constant bias of 222 /iV across the QPCs. The two data 
sets are generally in good agreement, with small system- 



atic deviations on the sides of the Coulomb peak and 
a certain scatter due to low statistics in the tails. The 
observed asymmetry between tunneling in and out pro- 
cesses (i.e., the difference in the maximum values of Fin 
and Tout) can be explained by the existence of a second 
degenerate quantum state in QD2. 

Having checked the consistency of the two methods in 
a regime where both are applicable, we test the correla- 
tion method in a regime with smaller signal levels. We do 
this by reducing the source-drain voltage on the QPCs. 
The step height A/''^' of the RTS is approximately pro- 

(n) 

portional to the bias whereas the noise level remains 
constant. The ratio of the two is the SNR relevant for the 
standard counting analysis, an insufficient SNR will re- 
sult in systematic measurement errors due to false counts, 
namely, an overestimation of the slower rate in case of an 
asymmetric RTS, or of both rates in case of a symmetric 
RTS. Assuming a certain current distribution of the am- 
plifier noise around the discrete current levels of the RTS, 
say a Gaussian distribution, the false count rate can be 
estimated as the number of statistically independent cur- 
rent measurements that lie outside a distance ctj AJ'^'^^/2 
from the mean. We can express it with the help of the er- 
ror function as 0.5[l-erf(SNR/2\/2)]/T„. The lower plot 
in Fig. m^a) shows a measurement of the signal-to-noise 
ratio along with the estimated false count rate calculated 
in this manner. The value for the SNR considered suffi- 
cient depends on the desired accuracy; here we require a 
SNR of more than 6 which results in a false count rate on 
the order of 10 Hz and which is reached for source-drain 
voltages larger than 150 /iV. 

In comparison, the measurement of Ftot shown in the 
upper plot of Fig. Ufa) demonstrates that the cross- 
correlation analysis is applicable down to significantly 
lower bias voltages, therefore reducing both the power 
dissipated by the sensors and the energy scale of the emit- 
ted radiation. As discussed, the best results are obtained 
close to the maximum of the peak where the rate is mea- 
sured reliably, i.e., with fluctuations below the statistical 
uncertainty due to the finite number of detected events, 
down to bias voltages of 22 /iV. Only below (and in the 
tails of the peak) the errors grow and eventually the anal- 
ysis algorithm fails. 

We now formulate a more precise criterion for compar- 
ing the two methods. In particular, it is first of all neces- 
sary to quantify the residual noise. For this purpose, we 
define as the standard deviation of the fluctuations 
in the function C(r) [cf. Fig. HJa)] measured in the ab- 
sence of a RTS signal. The ratio (cr^Vfi"^<''2"V^^ can 
be considered as a measure for the success in suppress- 
ing the noise by current cross correlation. However, the 
quantitative meaning of the noise level in the correlation 
case is different compared to the counting case. The ac- 
tual parameter of interest is the measurement uncertainty 
caused by this noise. Calculating it in the general case is 
a nontrivial task, on the one hand, because of the com- 
plexity of the analysis algorithm and, on the other hand. 
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FIG. 4: (Color online) (a) Top: tunneling rate Ftot between 
QD2 and lead as a function of QPC bias (applied to both PCI 
and PC2) and S gate voltage determined by cross-correlation 
analysis. Bottom: SNR for the counting analysis and ex- 
pected false count rate as a function of QPC bias. The 
false count rate was calculated from the SNR data assum- 
ing Gaussian amplifier noise with an autocorrelation time of 
T„ = 0.1ms. A SNR of 6 will result in a false count rate 
on the order of 10 Hz and can therefore be considered as the 
minimum requirement for the counting analysis. It is reached 
for bias voltages above approximately 150 ^V. (b) Examples 
of the time dependence (left column) and current distribution 
function P{I) (right column) of I2 recorded at three different 
QPC bias voltages indicated in the top graph in (a) . The RTS 
as a component of the current is recognized by the naked eye 
in all three cases, but only the trace "I" allows for a determi- 
nation of the transition rates without significant error when 
analyzed with a counting algorithm. 



because of the many experimental variables that play a 
role such as the absolute value of Ftot, R-TS asymmetry, 
measurement bandwidth, noise spectrum, and differences 
between the two channels (i.e., in the parameters and 

crj"' ) . We therefore restrict our discussion to the specific 
measurement situation discussed in this paper, in par- 
ticular, to the case of nearly identically coupled QPCs 
(ai = a2 ^ 1). We ask this question: by how much, 
starting from the limiting counting SNR of 6, can we 
reduce the signal strength A/ until we expect the corre- 
lation procedure to generate the same absolute error of 
about 10 Hz in Ftot? We write this "figure of merit" as 



AIrr 



AL 



min.X 



An) 



min.C 



J") 



(n) (n) A/min.C 



(8) 



The third factor in Eq. ([8]) is the original (inverse) SNR 
for the counting algorithm. The first factor can be consid- 
ered as the analog for the cross-correlation case, relating 
the signal strength to the residual noise cr^^ in C(r). It 
was determined with a numerical simulation. In apply- 
ing the data analysis algorithm to randomly generated 
time traces imitating the experimental ones (symmetric 
RTS with overlaid Gaussian noise, low-pass filtering with 
3 kHz, Ftot = 0.3 kHz, and T — 5s), the measurement un- 
certainty is obtained from the scatter in the output. The 
minimum AI for an error below 10 Hz determined in this 
way was given by 11 (u^''')^/^. Finally, the second factor 
in Eq. ^ is the noise reduction achieved in experiment; 



we measured cr^^ « 2.2 x lO'^^A^, a["' « 21 x IQ-^^ A, 



(n) 



and a. 
find 



(") 



16 X 10 A. Plugging in these numbers we 



AL 



min.X 



AL 
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i = 0.15. 



min.C 



(9) 



This means that in the case of the correlation experiment 
one can obtain meaningful values for the tunneling rates 
for signal-to-noise ratios approaching 1. 

To summarize, we have measured charge fluctuations 
on a GaAs DQD in a time-resolved manner simulta- 
neously with two QPC charge sensors. By evaluating 
their cross-correlation function and third-order correla- 
tors, we are able to determine the two time constants 
of tunneling back and forth between one dot and the 
adjacent lead. Obtaining the same information directly 
from either of the two QPC signals requires a signifi- 
cantly larger RTS amplitude because of the limitation 
due to amplifier noise. An interesting prospect is the ap- 
plication of the correlation technique to radio-frequency 
QPC setups2ii2Sii2^ where it would allow us to push the 
shot-noise limitation to the detection bandwidth toward 
the regime of charge qubit coherence times. 
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